The inequality of charge and spin diffusion coefficients 

Sandipan Prainanik and Supriyo Bandyopadhyay * 

'66' 

Q ' Department of Electrical and Computer Engineering, 

O 

CN I Virginia Commonwealth University, Richmond, VA 23284, USA 

Marc Cahay 

Department of Electrical and Computer Engineering and Computer Science, 

University of Cincinnati, Cincinnati, OH 45221, USA 

^ • Abstract 

1^ ■ Since spin and charge are both carried by electrons (or holes) in a solid, it is natural to assume 



C^ 



CJ 



that charge and spin diffusion coefficients will be the same. Drift-diffusion models of spin transport 
typically assume so. Here, we show analytically that the two diffusion coefficients can be vastly 



'^ ■ different in quantum wires. Although we do not consider quantum wells or bulk systems, it is 

O . likely that the two coefficients will be different in those systems as well. Thus, it is important to 



distinguish between them in transport models, particularly those applied to quantum wire based 



^ ■ devices. 
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In the drift-diffusion model of spin transport, it is customary to assume tliat tlie same 



diffusion coefficient 'D' describes charge and spin__diffusion. T 
in the hterature (see, for example, refs. [l|, |2 



iMU). Ref. 



lis assumption is commonplace 
6| considers a two dimensional 



system with different spin and charge diffusion coefficients but ultimately assumes that the 
bare spin diffusion coefficient is the same as the charge diffusion coefficient. Ref. |7| also 
examines this issue, and based on an heuristic assumption that spin transport is analogous 
to bipolar charge transport, reaches the conclusion that the two diffusion coefficients are 
equal as long as the populations of upspin and downspin carriers are equal. In spin polarized 
transport, the two populations are unequal by definition. Therefore, it is imperative to 
examine if these two diffusion coefficients are still equal in spin polarized transport, and 
if not, then how unequal they can be. In this paper, we show that these two diffusion 
coefficients can be vastly different in quantum wires. Although we do not consider quantum 
wells and bulk systems, there is no reason to believe apriori that even in those systems, the 
two diffusion coefficients will be equal. 

We first consider a narrow semiconductor quantum wire where only the lowest subband 
is occupied by carriers at all times. All higher subbands are unoccupied. We will assume 
that there are Rashba |8j] and Dresselhaus 91] spin orbit interactions in the wire, but no 
external magnetic field to cause spin mixing [1(|. In that case, we can ignore the Elliott- 
Yafet spin relaxation mechanism ll| since it will be very weak. Spin relaxation via hyperfine 
interaction with nuclear spins, or via the Bir-Aronov-Pikus mechanism 12|], is also typically 
very weak in semiconductor quantum wires with only one kind of carriers (electrons or holes, 
but not both). Therefore, the only spin relaxation mechanism that is important is the 
D'yakonov-Perel' relaxation |l3l |. 

In the single channeled quantum wire, we will prove two remarkable results: (i) spin will 
relax in time (i.e. the spin relaxation time r^ will be finite), but it will not relax in space (i.e. 
the spin relaxation length Lg will be infinite), and (ii) if the drift-diffusion model is valid in 



this system (this model relates Lg and r^ as Lg = y/D/Tg), then we must conclude that the 
spin diffusion coefficient Dg is infinite. However, since there is scattering in the system, the 
charge diffusion coefficient Dc must be finite. Therefore, the two diffusion coefficients are 
completely different. This is an extreme case, but even in less extreme cases (multi-channeled 
quantum wires), these two coefficients can be very different. Below, we provide an analytical 
proof for the single channeled quantum wire case. 



Consider an ensemble of electrons injected in a quantum wire at time t = from the end 
X = as shown in Fig. 1. Only the lowest subband is occupied in the wire at all times. There 
is an electric field E^ driving charge transport, and there is also a transverse electric field 
Ry breaking structural inversion symmetry, thereby causing a Rashba spin orbit interaction 
[Sl]. We will assume that the quantum wire axis is along the [100] crystallographic direction 
and that there is crystallographic inversion asymmetry along this direction giving rise to 
Dresselhaus spin-orbit interaction [9|. We choose this system because it is the simplest and 
includes the two major types of spin orbit interactions found in semiconductor nanostruc- 
tures, namely the Rashba and the Dresselhaus interactions. Ref. [5| has considered this 
system within the framework of the drift- diffusion model and shown that there is a single 
time constant describing spin relaxation. In contrast, spin relaxation in a two-dimensional 
system (quantum well) may be described by more than one time constant [Sj. 

For illustration purposes, we will assume hypothetically that the spin injection efficiency 
is 100%, so that at x = 0, all electrons are spin polarized along some particular, though 
arbitrary, direction f/o in space. Their injection velocities are not necessarily the same (in 
fact, they will be drawn from the Fermi-Dirac distribution in the contact). We are interested 
in finding out how the net spin polarization of the ensemble (|(S)|) decays in time or space 
due to the D'yakonov-Perel' process. 

In the quantum wire, the electrons experience various momentum relaxing scattering 
events. Between successive scattering events, they undergo free flight and during this time, 
their spins precess about a velocity-dependent pseudo-magnetic field Bso(fx) caused by 
Rashba and Dresselhaus spin-orbit interactions. This magnetic field can be shown to be 
spin-independent. 

The spin precession of every single electron occurs according to the well-known Larmor 
equation: 

§ = ^(v.) X S, (1) 

where S is the spin polarization vector of the electron and fl{vx) is a vector whose mag- 
nitude is the angular frequency of spin precession. It is related to Bso(^a;) as ft{vx) = 
{giJ,B/h)'Bso{vx), where g is the Lande g-factor in the material and /ib is the Bohr magne- 
ton. This is actually the well-known equation for Larmor spin precession and can be derived 
rigorously from the Ehrenfest Theorem of quantum mechanics. 



The vector fllv^) has two contributions due to Dresselhaus and Rashba interactions: 



Q(v:r) = fl-oivx) + 0-r(v^ 



iRlf^zJ 



(2) 



where the first term is the Dresselhaus and the second term is the Rashba contribution. 
These two contributions are given by 
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VxX = CdoVxX 



EyV^Z = CroVxZ, 



(3) 



where Wz, Wy are the transverse dimensions of the wire, 042 and a^^ are material constants, 
X is the unit vector along the x-direction and z is the unit vector along the z-direction. 

Note that the vector ft lies in the x-z plane and subtends an angle 6 with the ±x-axis 
(quantum wire axis) given by 



6 = arctan 
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(4) 



Note also that since 9 is independent of Vx-, the axis (but not the magnitude) of both 
Q, and Bgo is independent of electron velocity. Therefore, every electron, regardless of its 
velocity, precesses about the same axis, as long as only one subband is occupied. The direc- 
tion of precession (clockwise or counter-clockwise) depends on the sign of the velocity and 
therefore can change if the velocity changes sign, but the precession axis remains unchanged. 
However, the precession frequency depends on the velocity and is therefore different for dif- 
ferent electrons as long as there is a spread in their velocities caused by varying injection 
conditions or random scattering. As a result, at any given instant of time t = to, the spins 
of different electrons will be pointing in different directions because they have precessed by 
different angles since the initial injection. Consequently, when we ensemble average over all 
electrons, the quantity | (S) | decays in time, leading to spin relaxation in time. 

To show this more clearly, we start from Equation ([1]) describing the spin precession of 



any one arbitrary electron: 
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= -X [nR{v.^)Sy] - y [Qd{vx)Sz - Qr{vx)Sx] + z ^^{v^x)^, 

where S'„ is the spin component along the ra-axis of that arbitrary electron. 
Equating each Cartesian component separately, we get: 

—QmVx^yi 

= CrqVxSx — CdqVxSz, 

CoO'^xSy 
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(5) 
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dSz 

dt 



(6) 



If every electron in an ensemble had the same Vx at every instant of time (no dispersion 
in velocity), then the last equation tells us that every electron would have the exact same 
spin components Sx, Sy and Sz at any instant of time as long as they were all injected at 
time t = with the same spin polarization. In that case, we could replace 5*^ in the last 
equation by the ensemble averaged value < Sn > over the entire ensemble, so that 
d\<S>\^ d < Sx >^ . d < S„ >^ d < Sz >2 



dt 



d < Sx >'^ d < Sy>'^ 
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2<Sx>^^^ + 2<Sy>^^^-,2<Sz>^^^ 
dt dt dt 

-2(roVx < Sy >< Sx > +2CroVx < Sx >< Sy > -2(doVx < Sz >< Sy > 
+2Cdo'Vx < SyX Sz> 



0. 



(7) 



In that case, | < S* > | will not decay in time and there will be no D'yakonov-Perel' spin 
relaxation in time. However, if Vx is different for different electrons either due to different 
injection conditions, or because of scattering, then we cannot replace Sn with < 5*^ > in 
Equation ([H]). As a result. Equation ([7]) will not hold, so that d\ < S > \/dt ^ 0, and there 
will be a D'yakonov-Perel' relaxation in time. As a result, the spin relaxation time r^ will 
be finite. 



Next, let us consider D'yakonov-Perel' spin relaxation in space. From Equation ([H]), we 
obtain (using the chain rule of differentiation) 

ciox dx dox ^ ^ 

dx dt dx 

doy dx doy /- c /- a 

—, rr — ~i — Vx — (^rqVxOx — ^^dqVxOz, 

dx dt dx 

dbz dx dbz ^ „ ,^s 

The above equation shows that the spatial rates dSn/dx are independent oi velocity. This 
is a remarkable result with remarkable consequence. It tells us that even if different electrons 
have different velocities, as long as they were all injected with the same spin polarization 
at X = 0, they will all have the exact same spin polarization at any arbitrary location x 
= Xq\ That is, every electron's spin at x = Xq is pointing in exactly the same direction. 
Therefore, we can always replace S'„ in the above equation by its ensemble averaged value 
< Sn > whether or not there is scattering causing a spread in the electron velocity between 
different members of the ensemble. Consequently, 

d\<S>\^ _ d<Sx>^ d<Sy>^ d<Sz>^ 



\JjtXj \AjtXj \Aj<Xj VajiAj 

= 2<Sx>^^^ + 2<Sy>^^^ + 2<Sz>^^^ 

lAliAj (XiXy iAltAJ 

= -2C/J0 < Sy >< Sx > +2C/J0 < Sx >< Sy > -2Cdo < Sz >< Sy > 

+2Cdo < SyX Sz> 
= 0. (9) 

Thus, there is never any D'yakonov-Perel relaxation in space as long as a single subband is 
occupied. Therefore, the spin relaxation length Lg is infinite. This is true whether or not 
there is scattering. 

The above result has been confirmed independently with a many-particle Monte Carlo 



simulation of spin transport in a single channeled quantum wire IGj. Here, we have provided 
an analytical proof. 

The foregoing analysis also shows that in a quantum wire with single subband occupancy 
and D'yakonov-Perel' as the only spin relaxation mechanism, there is a fundamental differ- 
ence between spin relaxation in time and spin relaxation in space. Spin can relax in time 
while not relaxing in space. The physical origin of this difference is explained below: 



From Equation ([3]), we see that the precession frequency for any arbitrary electron is 
given by 

^ = Mit) = ^Jcho + O^it) = Cov.it), (10) 

where 0(t) is the angle by which the electron's spin precesses in time t. 

If all electrons are injected with the same spin polarization at time t = 0, then the angle 
by which any given electron's spin has precessed at time t = to is 

0(io) = Co / " v.{t)dt = Co[x(to) - x(0)] = (do, (11) 

where do is the distance between the location of the electron at time to and the point 
of injection. Obviously do is history- dependent, because different electrons with different 
injection velocities and/or scattering histories would traverse different distances in time to- 
Consequently, if we denote the angle by which the n-th electron's spin has precessed in time 
to as 0„(to), then 0i(to) ^ ^2(^0) ¥" •••0m(^o)- As a result, if we take a snapshot at to, we 
will find that the spin polarization vectors of different electrons are pointing in different 
directions. Therefore, ensemble averaged spin at to is less than what it was at time t = 0. 
Consequently, spin depolarizes with time leading to temporal D'yakonov-Perel' relaxation. 
The spatial rate of precession, on the other hand, is obtained as 

d(}){t) d(j){x) dx d(f){x) 
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Co- (12) 



Therefore, the angle by which any given electron's spin has precessed when it arrives at a 
location x = Xq is 

0(Xo) = p ^dx = Co P dx = CoXo. (13) 

This angle is obviously history-independent since it depends only on the coordinate Xq which 
is the same for all electrons at location Xq, regardless of how and when they arrived at that 
location. In fact, an electron may have visited the location Xq earlier, gone past it, and then 
scattered back to Xq. Or it may have arrived at Xq for the first time. It does not matter. 
The angle by which an electron's spin has precessed when it is located at Xq is a constant 
independent of past history. Therefore, if all electrons were injected with their spins exactly 
parallel to each other at a; = 0, then every single electron at x = Xq has its spin polarization 
vector pointing in the same direction as every other electron, and the ensemble averaged 



magnitude of spin at x = Xq is the same as that at x = 0. Consequently, spin does not 
depolarize in space and there is no D'yakonov-Perel spin relaxation in space, unlike time. 

Since spin relaxes in time but not in space, the relaxation time (r^) is finite whereas the 
relaxation length (Lg) is infinite. According to the drift-diffusion model, these two quantities 
are always related in steady state as [^ 

L, = ^/d^,, (14) 

where Dg is the spin diffusion coefficient. Note that the quantities Lg Dg and Tg are spin 
transport constants. As such, they are independent of both space and time. 

Since Lg is infinite while Tg is finite, the only way the above equation can be satisfied 
is if the steady-state spin diffusion coefficient Dg is infinite. But the steady state diffusion 
coefficient Dc associated with charge transport is certainly finite since we have frequent 
momentum relaxing scattering in our system. Therefore, there must be two very different 
diffusion coefficients Dg and Dc associated with spin and charge diffusion. This completes 
our analytical proof that Dc ^ Dg. 

Two final questions remain regarding the generality of the above result. First, is it only 
valid for the extreme case of a quantum wire with single subband occupancy (single channeled 
transport) and second, is it only true for Dyakonov-Perel relaxation? We cannot treat the 



case of multi-channeled transport analytical 
using Monte Carlo simulation in both space 



y, but we have examined that case numerically 



14j and time ISj . We studied spin transport in 



a GaAs quantum wire of cross section 30 nm x 4 nm, where multiple subbands are occupied 
and Dyakonov-Perel' relaxation does occur in both time and space. At a lattice temperature 
of 77 K and a driving electric field E^ = 2 kV/cm, the value of Lg extracted from that study 
is ~ 10 fim while the value of r^ ~ 1 nsec. This yields Dg ~ 10'^ cm^/s (from Equation flT^ ) 
which is still several orders of magnitude higher than the charge diffusion coefficient Dc in 
the same quantum wire calculated under the same conditions 171, [iSj. Thus Dg ^ D^ even 



in multi-channeled transport, and the two quantities can be vastly different. 

Finally, what if we include other modes of spin relaxation, such as Elliott-Yafet [1 ll] ? 
If Elliott-Yafet is the dominant mode, then spin relaxation is intimately connected with 
momentum relaxation. In that case, the charge diffusion constant, determined by momentum 
relaxing scattering, and spin diffusion constant may not be as unequal. Nontheless, there is 
no reason to assume apriori that the two diffusion coefficients are exactly equal even in this 



case. A rigorous Monte Carlo simulation (based on random walk model) recently carried 
out by us has shown that the two diffusion coefficients, in general, are vastly different. How 
different they are depends on the details of the scattering processes that relax momentum 



and spin 2G|. 



In conclusion, we have shown that in quantum wires, the spin and charge diffusion coeffi- 
cients are vastly different. Although we have not examined quantum wells and bulk systems 
in this study, there is no reason to pre-suppose that the charge and spin diffusion coefficients 
will be equal in these systems either. Thus, it is important to distinguish between these two 
diffusion coefficients in solid state systems. 
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Figure captions : 

Figure [H A quantum wire structure of witli rectangular cross section. A top gate (not 
drawn) applies a symmetry breaking electric field Ey to induce Rashba interaction. A battery 
(not drawn) applies an electric field —E^x, E^ > 0, along the channel. Spin polarized 
electrons are injected at x = 0. These electrons travel along x and may gradually lose their 
initial spin polarization. We investigate the spin depolarization of these electrons in time 
domain as well as in space domain. 



